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' We derive expression for energy flux in terms of lattice normal mode coordinates. Energy trans- 

mission across solid junctions from lattice dynamic point of view is given and its relation with 
atomic masses, lattice constants, and group velocities is clarified. A scattering boundary method 
(SBM) is proposed for calculating the amplitude transmission across solid junctions. The phonon 
transmission coefficients and thermal conductance are calculated for two kinds of acoustically mis- 
matched junctions: different chirality nanotubes (11,0) to (8,0), and Si-Ge interface structure. Our 
jy^ ' calculation shows a mode-dependent transmission in nanotube junction due to the high symmetry 

vibrating motions for nanotube atoms, indicating its possible important role in nanotube mixture 
thermal conductance. Energy transmission and Kapitza conductance across the Si-Ge interface [001] 
are calculated for the Si-Ge diamond-type structure. It is shown that the energy transmission across 
the Si-Ge interface depends on the incident angle and on the interface mode conversion. A critical 

. incident angle about 42° is numerically found for waves incident from Ge to Si. Our numerical 

result of the Kapitza conductance at temperature T = 200 K is Gk = 4.6 x 10 s WK -1 m -2 . We 
find numerically scaling law Gk oc T 2 ' 87 for [001] interface at low temperature. 



through the dynamic equations across the interface. 
Ref. has suggested calculating the phonon transport 
via Green's functions. However, the Green's function 
method gives only a formal solution for the scattering 
problem. A practical calculation procedure is not pro- 
vided. Ref. |8| used phonon mode of wave incident from 
one lead to calculate the energy transmission through 
the amplitudes of the transmitted/reflected waves by 
ikj = PAj2q v qzAl/(p B v 0z A%) without any justification. 
How to calculate the energy transmission from the am- 
plitudes of the transmitted/reflected waves, and what 
is its relationship with atom mass, lattice constant and 
group velocity? It's not clear. On the other hand, the 
well-known Landauer formula for electronic conduction 10 
can be generalized to thermal transport for such calcu- 
lations. Landauer formula for ballistic heat transport 
across junctions has been usedii*^ for the prediction of 
universal quantum heat conductance at very low temper- 
atures. The energy transmission formula derived under 
continuum assumption may not be applied straightfor- 
wardly to systems on nanoscale where atomic details are 
important. 

In this paper, we report our results of the energy trans- 
mission across solid junctions taking into account atomic 
details from lattice dynamic point of view both theoreti- 
cal and numerically. In Section II, we will first derive the 
energy flux by lattice dynamics method and then give 
the energy transmission across the interface involving 
atomic mass, lattice constant and group velocity, start- 
ing from the energy conservation law. An analytic result 
for ID acoustic mismatched chain model from our for- 
mula is also illustrated. We further propose computing 
the amplitude transmission coefficients by solving a set of 
dynamical equations with a scattering boundary condi- 
tion method (SBM) to cope with the complex interface. 
In Section III, we first compute the transmission coeffi- 



o 
o 



o 
o 



- 1—^ 

X 



PACS numbers: 66.70.+f, 44.10.-K, 05.45.-a 



I. INTRODUCTION 



Rapid progress in the synthesis and processing of ma- 
terials with structures of nanometer length scales has cre- 
ated a demand for the understanding of thermal trans- 
port in nano-scale low dimensional devices 1 - 2 - 3 - 4 . Recent 
experimental and theoretical studies have revealed novel 
features of phonon transport in these systems, such as 
the size-dependent anomalous heat conduction in one- 
dimensional (ID) chainS and the universal quantum ther- 
mal conductance^. Thermal transport in nanostructures 
may differ from the predictions of Fourier's law based on 
Q\ . bulk materials; this may happen because of the existence 
of many acoustically mismatched interfaces in nanostruc- 
tures and because the phonon mean free path is compa- 
rable to the size of the structure 4 . An understanding of 
the thermal conduction across acoustically mismatched 
solid interfaces is a necessary requirement for thermal 
Ph I transport engineering. 

The study of thermal transport across interfaces dates 
back as early as to 1940s when Kapitza resistance^ 
was reported and much work has been done in this 
field 4 ^. In general, theoretical modeling of this prob- 
lem has been undertaken either by the acoustic-mismatch 
model (AM) with scalar elastic waves, or by the 
diffuse-mismatch model (DM) with Boltzmann transport 
equation^. Some numerical methods such as molecular 
dynamic simulation 7 have also been used. While AM 
and DM models provide some useful reference calcula- 
tions, scalar wave model and Boltzmann transport equa- 
tion are only phcnomenological descriptions and they 
have ignored the complexity of the interface. Atomic- 
level lattice dynamic (LD) approach should be the right 
way of capturing the mechanism of heat transport. Au- 
thors in Ref. have proposed to calculate the ampli- 
tude transmission by connecting the neighboring atoms 
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cients of a Carbon-nanotube junction showing a mode- 
dependent transmission in nanotubc junction due to the 
high symmetry of vibration in nanotubes. Then for the Si 
and Gc diamond-typed structure by linearized empirical 
Tersoff potential, we calculate the energy transmission 
across the Si-Ge interface [001] to show its dependence 
on the incident angle. A critical incident angle about 
42° is numerically found for waves incident from Ge to 
Si. We also calculate the Kapitza conductance for Si- 
Ge interface at different temperatures and find its value 
G K = 4.6 x 10 s WK _1 m- 2 when T = 200 K and numer- 
ically get its scaling law T 2 87 at low temperature. 



II. THEORY 

A. Energy Flux 

In this section we derive the energy flux in terms of 
the normal mode coordinates of the atomic vibration 
from the energy continuity equation. We consider sys- 
tems with perfect leads on the left and right with an 
arbitrary interaction at the junction. The Hamiltonian 
for each lead system (as well as the whole system) of vi- 
brating atoms under linear approximation takes the form 



H = 



where I or I' denotes a unit cell, i,j the position in a 
cell, a, /3 the direction of vibrating motions of atoms, 
and uf i the displacement from equilibrium of the atom 
(I, i) with equilibrium position R;,j. The dynamic matrix 
eigenvalue problem can be written as 
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where D 



— Kff,,^*^'-^. The eigen- 

vector for Eq. (0 is ef„(q), where i denotes the posi- 
tion in a cell, a the direction of vibrating motion, n 
the branch of polarization. These eigenvectors satisfy 
Y,i, a e"„(q)e*"/(q) = $„,„'• So the vibrating motions 
of atoms can be expanded into the summation of these 
eigenvectors 
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where is a normal mode coordinate of the vi- 

bration, and N is the number of the unit cells. Each 
normal mode represents a single harmonic oscillation 
Qq(t) + w 2 (q)Qq(t) = 0. A local energy density can 
be defined through the energy in cell I from Eq. (JIJ, as 
p(r) = J2i £ ifi( r — Ri) where R; is a lattice vector. 

Next, an expression for the energy current in the z di- 
rection can be derived from the energy continuity equa- 
tion, dp/dt + V ■ j = 0, where j is the energy current 



density. We transform the energy continuity equation 
into the momentum space, 



(n(q,t)) t + ig*(J*(q,t)} t = 0, 
by the following Fourier transform 
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In these equations, f2(q, t) and J(q, t) are the energy den- 
sity and energy current density in momentum space, re- 
spectively, and the bracket ( ) t denotes the time average. 
Here we have assumed that the energy current propa- 
gates along z direction and so the time averaged energy 
current components along other two directions equal to 
zero. A relation of the total energy current in real space 
and in momentum space is 



Km J(q,t) = / j(r,t)d 3 r. 
J 



Using Eq. J3J and Eq. Q , we can get the average energy 
flux along the z direction, 
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where V is the volume of the lead part. 

With the help of Eq. Q, Eq. ©, and Eq. ©, the 
energy density in momentum space is given by 

tt(q ,t) = ]>> e -^- R ' = i Yl 

I q,q', n,n' 



i.a. 



(q)OWW,«}. (9) 



where Qq's are normal mode coordinates for the n-th 
branch phonons. Combining Eq. (JSJ and Eq. @, the 
time-averaged energy flux along z-direction is 
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Note that during the derivation of Eq. i|10|) . il(q, t) 
should be first differentiated with respect to time and 
then takes the limit q z — > 0. Eq. 110|) is the formula for 
the energy flux in terms of normal mode coordinates from 
lattice dynamic point of view. It provides the basic start- 
ing point for studying the energy transmission across the 
interface. 



B. Energy Transmission 

Once the energy flux formula Eq. (|10fl is obtained, we 
can further discuss the problem of energy transmission 
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across the interface. The idea for lattice dynamic ap- 
proach to study energy transmission across the interface 
was first proposed in Ref. 0- Assuming an eigen mode 
lattice wave for Eq. J2J incident from one lead, it will 
be refracted by the interface. The reflected and trans- 
mitted waves are assumed the linear combination of the 
eigen mode lattice waves of Eq. © at each side of the 
lead. These amplitude transmissions are solved through 
the equations of motion for the boundary atoms. But 
it should be noted that these amplitude transmissions 
are not energy transmissions. What are their relation 
to the energy transmissions? As already cited in the in- 
troduction part of this paper, several assumptions^ are 
made to calculate the energy transmissions from these 
amplitude transmissions. Apart from this problem, any 
eigenvector for Eq. (PJ can be multiplied by a whatever 
coefficient and it is still the eigenvector. It means that 
the choice of eigenvector is not unique because of the lin- 
ear nature of the harmonic system. This indeterminacy 
will complicate the calculation of the amplitude transmis- 
sions and makes the amplitude transmission coefficients 
non-unique. However the energy transmissions are deter- 
mined by the system and have unique values. The key to 
this problem is energy conservation, which means that 
the energy flux from the left and right lead should be 
equal. In this part, we use the energy flux Eq. 1)10(1 de- 
rived from the energy continuity equation in the previous 
section to clarify these issues. 

We assume that the vibration of the atoms for each 
incident/tranmitted normal mode is 



(11) 



where ef n (q) is the eigenvector for dynamic matrix 
Eq. © and it satisfies J2 ia e"„(q)e*"/ (q) = 5 n , n >. It 
is obvious that our choice for the eigen mode is unique 
and deterministic. Compared with the choice of incident 
normal mode in Ref. 0, the difference is that the atomic 
mass should be included. The reason for this is simple: 
the eigenvectors of Eq. J2J) are not the travelling wave so- 
lution for the system and its relation to travelling wave is 
to be divide by y / m7 according to the definition of Eq. J5J. 
We assume the same form, Eq. (|11|) . for all the incident, 
reflected, and transmitted waves. When one particular 
mode ufl (to, q) is incident from the left lead, it will 
be scattered by the interface. So the reflected wave is 
u), — q') and the transmitted wave 



given as \^ n , t„,„u 
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In these equations, t°,° is the 



amplitude transmission/reflection coefficients from mode 
n in the lead a = L to mode n' in lead a 1 = R, L. Note 
that q' and q" satisfies ui — ui n >(—q') — ui n >(q"). Similar 
expression can be written down for the right lead. The 
total wave for one particular lead, say the left lead, is an 
arbitrary superposition of all these eigen modes. Thus, 



the motion of the atoms is described by the wave: 
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where a a n is an arbitrary amplitude of the mode n in 
lead a. With the help of our newly derived energy flux 
Eq. I|10[l for lattice dynamics, after transforming Eq. 1|12|) 
into normal mode coordinates, we can get the heat flux 
for the lead a = L,R: 
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Since the energy must be conserved, there is no net en- 
ergy accumulation in the junction, which means that 
the time-averaged energy currents from both sides are 
equal, I L = I R , where the energy currents I a — I a ■ S a , 
a = L, R and S a represents the cross area for each lead. 
This condition leads to the following identity for the 
transmission amplitudes, 
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where f z is the length of unit cell along the z direction, 
v a n = duj°/dq z . This equation is analogous to the unitar- 
ity condition for the scattering matrix in electronic trans- 
mission. To clearly understand the meaning of Eq. H14fl . 
we can write it into a matrix product equation t^vt = v 
with 
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The above expression and the technique we used are sim- 
ilar to the continuum case^ but the mathematics in- 
volved in the lattice dynamic framework is simpler and 
clearer. Several interesting conclusions can be made from 
Eq. CU) and Eq. fTB) . 

(1) The amplitude transmissions from each lead are 
not independent. They must obey Eq. (|14|l so that the 
energy will be kept conserved. It means that which lead 
is chosen for calculating the energy transmission does not 
matter. They will give the same energy flux from cither 
of the lead. In practice, this should be a good criterion 



to check the validity of the algorithm in calculating the 
amplitude transmission. 

(2) From the diagonal terms of Eq. (|14[l. the energy 
transmission T°,° from mode (a, n), which means a lead 
and n branch, to the mode (cr',n') is given by 
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FIG. 1: Schematic show for scattering boundary condition. 



The total reflected and transmitted energy transmission 
for (cr, n) is given by 



thermal average, we can get the Kapitza conductance as 
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(3) The energy transmission's relation to atomic mass, 
lattice constant and the group velocity is explicitly ex- 
pressed in Eq. (|17af> . It can be seen that the atom 
mass rrii for each lattice site R| { should be considered 
in each eigen mode wave as in Eq. ifTTjl. Ref. used 
tkj = PAj2 q v qzAl/(pBVozAl) to calculate the energy 
transmission from the amplitude transmission A q . This 
equation is correct only if the atomic masses are the same 
in a unit cell. In actuality, the ratio of Silicon mass den- 
sity to Germanium mass density is about 0.4375 at tem- 
perature T = 300 K, while the accurate coefficient should 
be 28/72 w 0.389 if the energy is conserved and lattice 
constants are same for both leads. Ref. Il3l improved this 
by using the atomic mass ratio in the final energy trans- 
mission to keep energy conserved. But when the unit cell 
of system is composite, the choice of incident mode wave 
of Ref. 0is wrong. It should be chosen as our Eq. (JTTJ. So 
the results of thermal conductance between metals and 
BaF2 in Ref. are arguable using the method in Ref. 8. 
The lattice constant is also considered in Eq. i|17a|) , which 
may have significant contribution when the difference in 
the lattice constants for the left and right lead crystal is 
large. 

(4) Comparing with the continuum caseii, we see an 
extra factor of the lattice dimension in z direction, f z . 
This factor can be included in the amplitude transmis- 
sion if we redefine our normal mode as that in Eq. I|ll|) 
multiplied by v^cciij where V ce \\ is the volume of a unit 
cell. 

Next, we quantize Eq. (|10|) with 
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We get the quantization of the heat current I — 
FEq^^W^lllalfl^fl' Similarly after quantizing 
Eq. (|13fl with the help of the relation Eq. i|14|) , and taking 
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where T n (q,o; n ) = J2 n , \t% L n \ 2 v%/v% and f(w n ,T) is the 
Bose-Einstein distribution at temperature T. The dis- 
persion relation and group velocity refer to the left lead. 
In deriving Eq. I|19(l , we have used the fact that phonons 
in the lead obey the Bose-Einstein distribution. 



C. Scattering Boundary Method 

The central issue now is to have an efficient method to 
calculate the amplitude transmission coefficients across 
junctions at the atomistic level. Ref. has proposed 
solving the amplitude transmissions through connecting 
boundary atoms for abrupt solid interface. But this ap- 
pears to be a difficult problem for general complex inter- 
faces. As illustrated in Fig. ^ the central part of atoms 
is not known so the method proposed in Ref. will not 
work. Transfer matrix method can be used for simple 
ID models with a few atoms in the central part to obtain 
analytical solutions. But this method is not numerically 
stable for a long ID system (such as a disordered ID 
system) or 3D systems because the errors get amplified 
by continuous multiplications of a matrix with eigenval- 
ues larger than 1. In fact, we have tried transfer matrix 
method for the nanotube junction and found large nu- 
merical errors of 30% or more. 

Here we propose what we call the scattering bound- 
ary equation method (SBM) as a solution. Each atom 
in the system satisfies the dynamic equation, — miUj 2 Ui + 



U; 



0, where the 3x3 matrix Ku is 



the force constants between atom i and j. The bound- 
ary conditions are of the form u qin > +5Z n r„U- 9 .„ for the 
incoming and reflected waves and ^n u q' n f° r ^ ne ou ^" 
going waves, where u qin and u' q , n are the eigen modes on 
the left and right leads, while the reflection and transmis- 
sion coefficients r n and t n are unknown. These equations 
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FIG. 2: Schematic show for ID acoustically mismatched 
chain model. 



in matrix form are illustrated as 



from the right lead as 
1 
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(22b) 



As illustrated in the schematic Fig. [21 there are two 
atoms, labelled 2 and 3 in the figure, in the connect- 
ing part. The equations of motion for the two atoms are 
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where Ka = — m^ 2 + J2j Kf. These equations can be 
both analytically and numerically solved by the conven- 
tional method if the matrix is square. In some cases 
for higher dimensions, however, the boundary conditions 
are complicated and the number of equations may be 
larger than that of variables. But these equations are 
not linearly independent. Nevertheless, these equations 
are consistent and simultaneous. Under this condition 
the scattering boundary equations can still be solved nu- 
merically by the standard singular value decomposition 
method^. 



D. Analytic Results on ID Model 

To demonstrate the formulas obtained in the previous 
part and the scattering boundary equation method, we 
give some analytic results on a toy ID model. More real- 
istic material solid junctions, e.g. nanotube conjunction 
and Si-Ge interface, will be discussed in the next sec- 
tion. We consider two acoustically mismatched chains of 
different atomic masses mi,m2, connected with springs 
of different stiffness K%, K 2l and lattice constant a, b, on 
each side as illustrated in Fig. El In the central part they 
are connected by a spring with stiffness Kf. The dynamic 
matrix for the ID chain gives the eigenvector e; = e lqxi 
and w| = — cos(qa)). From Eq. the incident, 

reflected and transmitted waves from the left lead can be 
written as 




where, following the notation of Eq. I|12fl . t LL and t RL 
are the reflected and transmitted amplitudes. Similar 
expressions can be written down for the wave incident 



{- mi Lu 2 + Ki + K f )u 2 - K x ui- K f u 3 = 0,(23a) 
{-m 2 uj 2 + K f + K 2 )u 3 - K f u 2 - K 2 u A = 0.(23b) 

For wave incident from the left lead, the SBM equations 
can be expressed as 
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where K 22 — —miuj 2 +Ki+Kf, K 33 — —m 2 uj 2 +Kf+K 2 
and qi 1 q 2 satisfy that ujl(qi) = ur(Q2)- For simplicity, 
we only consider the condition that the two spring chains 
are connected directly by the spring Kf. If there are 
more atoms at the central part, we only need to add 
equations like Eq. I|23|l for these atoms to Eq. I|24|). But 
the scattering boundary condition in Eq. I|24|l . last four 
rows, does not change. From Eq. (|24|l . we can get 
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Similarly for wave incident from the right lead, we have 
SBM equations 
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From Eq. (|26|l . we can also get t RR and t LR . Next we 
construct Eq. (|15fl for ID model as 



t LL t LR 
h RL 4-RR 




(27) 



Matrix Eq. (|15fl . t^vt = v, gives four equations for this 
ID model, 
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These equations can be easily verified using the results 
from Eqs. (PH 1251 l2*S|l . These four equations also guar- 
antee the same energy transmission both for the left and 
the right lead, that is 
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T 1 
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It can be easily seen from the simple ID case that am- 
plitude transmissions from the left and right lead are not 
independent and should satisfy the relation (|14|l to keep 
the energy conserved. 

Furthermore, for ID quantum thermal energy cur- 
rent, with the help of the relation l|14|) . after quantizing 
Eq. we get the formula 
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where T n (uj) — \ tmn\ ^n/^m i s tri e energy transmis- 

m 

sion probability, f(uj,T a ) is the Bose- Einstein distribu- 
tion for the left or right lead. Eq. (|30() is the Landauer 
formula for quantum thermal energy flow in one dimen- 
sion. 



III. SIMULATION 

In this section, we report our results for two kinds 
of solid junctions: nanotube junction and the Si-Ge 
interface. Nanotube has been found to be good can- 
didate for thermal conduction materials and the ther- 
mal conduction has been improved through adding nan- 
otube mixture^. But the interface between nanotubes 
are inevitable. How they will affect thermal conduction 
becomes an interesting problem. Thermal conduction 
across the Si-Ge interface has been studied by Ref. 
using the simplified fee lattice model and recently by 
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FIG. 3: (a) Structure of a (11,0) and (8,0) nanotube junction, 
(b) The energy transition T n (oj) as a function of angular fre- 
quency Lu. 

Ref. El in diamond structure with equal lattice constants 
for Si and Ge and with fee neighboring atoms considered. 
In our paper, we will consider more realistic crystal struc- 
ture for Si and Ge to investigate the energy transmission 
dependence on incident angle and low-temperature scal- 
ing behavior of Kapitza conductance. 



A. Nanotube Junction 

The semiconductor nanotube junction structure (11,0) 
and (8,0) is first constructed by a geometrical method as 
in Ref. 16. However there are many defects in the geomet- 
rically constructed structure. Next we optimize the junc- 
tion structure by a second-generation Brenner potential 1 ^ 
to let the atoms get their equilibrium positions. The op- 
timized structure is shown in Fig.^a). Then under small 
displacement, the force constants are derived numeri- 
cally from the same potential. The phonon dispersion 
calculated from these linearized force constants for nan- 
otube (11,0) is illustrated in Fig.0J in which four acoustic 
branches are considered for energy transport: the longi- 
tudinal mode (LA), doubly degenerate transverse mode 
(TA), and the unique twist mode (TW) in nanotubes. 

Following the SBM method, the SBM equations are 
constructed and it involves 504 equations and 455 vari- 
ables. However, some equations are not linearly indepen- 
dent and are consistent because the motion for each atom 
is connected with the motions of others forming travel- 
ling waves and they are confined together to construct 
the nanotube structure. We also numerically verified the 
consistency of these equations through the calculation 
of ranks for SBM equations. The rank we get is equal 
to 455 within numerical accuracy, which is equal to the 
number of variables. These equations are solved through 
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q T/rc 

FIG. 4: Phonon dispersion for nanotube (11,0). The wave 
number q is expressed in terms of the length of nanotube 
translation vector T. 



the standard singular value decomposition method^*. 

We then discuss the SBM results of the transmission 
coefficients for a nanotube junction, shown in Fig.[!|Jb). 
Although all modes of a given frequency are considered, 
we did not find mode-mixing or mode-conversion behav- 
ior among acoustic modes at the lower frequency range. 
For example, the reflected and the transmitted waves 
across the junction for the incident LA mode waves are 
both only LA modes. We think that this is due to the 
high symmetrical properties of atomic motion for nan- 
otubes. Each travelling wave on nanotube has its own 
symmetrical property resulting from the symmetry of the 
nanotube structure. The symmetry on the left and right 
lead is different. It is hard to convert one specific high 
symmetrical motion to the other different symmetrical 
motion through the conjunction except when the sym- 
metrical properties on each lead are the same. We think 
that this kind of mismatch in the symmetry of motion 
for nanotubes will play an important role in the thermal 
conduction of nanotube mixtures. 

The LA modes are common symmetrical motion for 
both the left and right lead. The transmission for LA 
mode stays around 0.8 with only small changes. This 
value is below the AM model prediction^ of 0.98 with the 
longitudinal group velocity 20.18 kms -1 and 20.95 kms -1 
for (11,0) and (8,0). In contrast, the transmissions of the 
TW mode and many other optical modes are nearly zero 
or very small. This appears related to the difference in 
rotational symmetries of these modes. The transmission 
for TA mode decreases with frequency. This can be ac- 
counted for two reasons: (1) Although the TA mode is a 
common symmetry for both left and right lead, the kind 
of transverse symmetry is destroyed by the central junc- 
tion part. It can be seen from Fig. [3] that the central 
junction part does not have transverse symmetry. (2) 
There is a mismatch of dispersion relation for TA mode 
because the TA mode has a nearly quadratic dispersion 



relation, as illustrated in Fig.|2] These reasons will affect 
the transmission for TA mode at high frequency when 
the vibration motion of atoms is fast. 

So it can be seen from above that energy transmis- 
sion across nanotube junctions is strongly dependent on 
symmetry property of the lead nanotube. This symme- 
try dependence is also observed by molecular dynamics. 19 
Usually such mismatch in symmetry will reduce the num- 
ber of possible transmission modes and thus blocks the 
thermal conduction across the junction. We also propose 
that this kind of mode-dependent transmission behavior 
may be important for further application such as phonon 
filters. 



B. Si-Ge Interface 

The Kapitza conductance across the Si-Ge interface 
have been calculated by lattice dynamic method in Ref. 8 
usingthe simplified fee lattice model and recently by 
Ref. [l3| in diamond structure with fee neighboring atoms 
considered. However, the relation of energy transmission 
to the incident wave angle and the temperature depen- 
dence of Kapitza conductance at lower temperature re- 
mains unclear. Here we use our newly derived formulas 
combined with the simulation to get an understanding of 
these problems in more realistic Si-Ge structure. We first 
discuss the energy transmission's relation to the incident 
angles and to wave mode conversion at the interface. The 
Kapitza conductance across Si and Ge interface is also 
calculated and its scaling law in relation to temperature 
at low-temperature regime is numerically analyzed. 

The lattice structure for Si and Ge crystal we con- 
sidered is the diamond structure. The atomic mass for 
Si and Ge are 28 amu and 72.61 amu. The lattice con- 
stants for each conventional unit cell are a — 5.43 A and 
a = 5.658 A, respectively. We use the Tersoff potential^ 
under small displacement to get the linearized force con- 
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FIG. 5: Phonon dispersion along T — L for Si and Ge. The 
thick lines are for Si and the thin lines for Ge. The wave num- 
ber is in terms of the maximum value along r — L direction. 
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stants. By Tersoff potential truncation function^ only 
the nearest 4 atoms are considered for each atom with 
the distance y/3a/4. We choose the primary cell consist- 
ing of two atoms (0, 0, 0) and (a/4, a/4, a/4) to construct 
the dynamic matrix. The phonon dispersion for Si and 
Ge along the T — L direction can be calculated through 
the linearized force constants from the dynamic matrix 
and are illustrated in Fig. [S] The maximum LA branch 
along r — L frequency for Si and Ge are 343 cm -1 and 
196 cm -1 . It can be seen that the linearized force con- 
stants for the nearest atoms well reproduce the Si and 
Ge phonon dispersion^ The group velocity is calculated 
through 



1 etfe 



2lo et 



(31) 



where D is the dynamic matrix and e the eigenvector 
as shown in Eq. (J2J). The group velocity of LA branch 
phonons along the F — L direction for Si and Ge is about 
7.8Km/s and 4.4Km/s. 

For the Si or Ge diamond structure, there are six modes 
of vibration as illustrated in Fig. [5] , among which three 
branches are acoustic and other three ones are optical. 
Using Eq. when one mode of wave u" in (uj,q) = 

-^=ef rl (q)e l< - q Ri ~" t ' ) is incident from the lead, it has 
three modes of reflected wave and three modes of trans- 
mitted wave. Due to linearity of the system, the fre- 
quency does not change during refraction. We choose 
the Si-Ge interface with normal in [001] direction. Since 
the system is homogeneous in the x and y direction, the 
reflected waves and the transmitted waves have the same 
momentum in these directions q' xl q" — q x \ q' y , q'y — q y . q' z 
and q" is found to satisfy w(q) = wsi(q') == w Gc(q")- K 
q' z or g" is complex, Im[g^] or Im[g"] should be negative. 
The force constants for the Si-Ge interface are chosen 
as Kij — (Kf* + K^ e )/2. Each atom in the incident 
lead part has the motion u re fl ec t = u(q) + riii(— q'jj + 
?~2u(— q 2 ) + r^u(~ q^); similarly for the motion in the 
transmitted wave Utransmit = iiu(qi')+<2u(q 2 ')+<3u(q3)- 
Connecting the Si-Ge interface provides us two equations 
of motion, 

( m Si u 2 • ^K,,)u, )jK,,u, = 0, (32) 

3 3 

(-m Ge w 2 +^K !:) -)u,;-^K IJ Uj = 0. (33) 

3 3 

There are all 8 different atoms in the above two equa- 
tions. Each atom should have the form of solution of 
the reflected u rc fl ct or the transmitted u transm it, which 
depends on its position. The scattering boundary equa- 
tions are constructed. There are 30 equations (two equa- 
tions of motion give 6 and another 8 atoms taking the 
scattering boundary condition solution give 24) and 30 
variables (8 atoms give 24 and 6 for reflected coefficients 
r and transmitted coefficients t) in the scattering bound- 
ary equations. This set of equations can be solved by a 



conventional method. When r and t is obtained, the en- 
ergy reflection coefficients R and transmission coefficients 
T for each mode are given by summation of Eq. Ijl6|l . 

We first report the result of the dependence of energy 
transmission on the incident angle. A continuum wave 
incident on the surface whose wave vector makes an angle 
9 with the surface normal is refracted in accordance with 
Snell's law. Does the refraction of phonon incident on 
the solid interface observe the Snell's law? Is there any 
new character for phonon refraction? We calculated the 
energy transmission for LA mode incident from Si to Ge 
and incident from Ge to Si. The results are illustrated in 
Fig. © and Fig.d 

It can be seen from Fig.Elthat when the angle for the 
phonon incident from Si to Ge increases, energy trans- 
mission decreases slowly. But once the angle is larger 
than 80°, the transmission decreases rapidly. In con- 
trast to the transmission from Si to Ge, LA mode wave 
transmission from Ge to Si shows a rich character as il- 
lustrated in Fig. For wave with the incident angle 
0^0, the transmission first decreases with the increase 
of frequency. But when the frequency goes over a cer- 
tain value, for example u> — 93 cm -1 for 8 = 36°, the 
transmission begins to increase. This behavior can be 
understood by the mode conversion at the interface. It 
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FIG. 6: Dependence of energy transmission on the incident 
angle. LA phonon incident from Si to Ge. 
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can be seen from Fig. [5] that the maximum frequency for 
TA modes in Ge is about 98 cm -1 . When the frequency 
is below this value, there are TA modes for the reflected 
wave; but over this value, the reflected wave can not be 
converted into TA modes. So due to lack of reflected 
modes, the transmission increases. 

Furthermore, for LA waves incident from Ge to Si, 
there exists a critical value, about 42°, above which 
transmission equals 0. We can estimate the critical an- 
gle with the help of SnelPs law for continuum wave. 
The group velocities for normal incident waves in direc- 
tion [001] are m 6.87 Km/s and vf e m 3.78 Km/s. 
The SnelPs law gives the critical angle from Ge to Si 
6 C = Bhr 1 (vf e /v% i ) = 33.4°. This value is a little lower 
than the observed critical angle value. Nevertheless, we 
think that the Snell's law gives a rough estimation of the 
critical value for the critical incident angle for phonons. 

When the incident waves from Si to Ge and from Ge 
to Si are both normal to the surface (9 — 0), the en- 
ergy transmission from both sides are about the same 
0.98. We use the acoustic mismatch modeti& to estimate 
4:ZsiZ Ge /(Z S i + Z Ge ) 2 w 0.97. It can be seen that the 
acoustic mismatch model describes the transmission for 
normal incident waves well. 

The Kapitza conductance with change of temperature 
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FIG. 7: Dependence of energy transmission on the incident 
angle. LA phonon incident from Ge to Si. 



is calculated using Eq. I|19|) . The results for Kapitza 
conductance and heat capacity using our linearized force 
constants are illustrated in Fig. |Hl The Kapitza conduc- 
tance for Si-Ge [100] interface calculated from our model 
is G K = 4.6 x 10 8 WK^m- 2 when T = 200 K. When the 
temperature goes over 200 K, we find that the Kapitza 
conductance changes little with the temperature and is 
saturated. For comparison, we plotted the heat capac- 
ity in Fig. [SJ It can be seen from the figure that the 
heat capacity continues to increase with the temperature 
when T > 200 K. In comparison with the heat capacity, 
the saturation of Kapitza conductance can be accounted 
by little contribution of energy transmission from high 
frequencies. At low temperatures, as illustrated in in- 
set of Fig. HO the Kapitza conductance scales as T 2 - 87 , 
while the heat capacity scales as T 3 in accordance with 
the Debye model. We have sampled enough points in the 
first Brillouin zone to ensure that Kapitza conductance 
and heat capacity converge numerically. However, due 
to the small deviation from 3, we cannot rule out the 
possibility that the exponent for Kapitza conduction is 
also 3. Ref. reported that Kapitza conductance scales 
as T 3 at low temperatures for fee interface irrespective 
of the properties for the left and right lead. The Kapitza 
conductance's dependence on temperature is an intrigu- 
ing problem^, with much experimental work done in this 
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field. Most experiments gave T a with a < 3 for solid 
interface as reviewed in Ref. |5|- So far no experimental 
result is available for the temperature dependence of the 
Si-Ge interface. What value should it be is an interest- 
ing problem. Compared with the results of Ref. |8|, our 
discrepancy from T 3 comes from the anisotropy of the 
energy transmission because of the diamond structure is 
used for the calculation of the transmission, while Ref. 
took the isotropic assumption in their calculation. 

IV. CONCLUSION 

In summary, we derive expressions for energy flux in 
terms of lattice normal mode coordinates. Energy trans- 
mission across solid junctions from lattice dynamic point 
of view is given and its relation with atom masses, lat- 
tice constants, and group velocities is also clarified. A 
scattering boundary method (SBM) is proposed for cal- 
culating the amplitude transmission across solid junction. 
Our calculation shows a mode-dependent transmission in 
nanotube junction due to the high symmetry vibrating 
motion for nanotube atoms, indicating its possible im- 
portant role in nanotube mixture thermal conductance. 
Energy transmission and Kapitza conductance across the 
Si-Ge interface are calculated for the Si-Ge diamond- 
type structure from linearized Tersoff potential. It is 
shown that the energy transmission across the Si-Ge in- 
terface depends on the incident angle and on the inter- 
face mode conversion. A critical incident angle about 



42° is numerically found for waves incident from Ge to 
Si. The critical angle in the reverse direction is much 
larger. The Kapitza conductance saturates to about 
G K = 4.6 x 10 8 WK _1 m- 2 for T > 200 K. We numer- 
ically get its scaling law T 2 87 for [001] interface at low 
temperature. 

In additions, we remark that 

(1) Eq. ijTil) and Eq. (|T5|l are the basic relations for 
amplitude transmissions in linear systems, irrespective of 
the details of junction part. The results for any algorithm 
to calculate amplitude transmission from lattice point of 
view should satisfy these relations to keep the energy 
conserved. It dose not matter which lead to be chosen 
for the energy transmission because the energy flux will 
be the same according to Eq. Ijl4|l . So Eq. (|14|l provides 
a criterion for checking the validity of the algorithm. 

(2) Lattice dynamic approach to solid junction is only 
under linear approximation, which may be accurate at 
low temperatures. However, for high temperature the 
nonlinear effect or anharmonic effect has to be consid- 
ered. How to include the effect of nonlinearity into 
phonon scattering at the solid junction will be an in- 
teresting problem. 
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